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Abstract 

The paper assesses stationary probability distributions in out-of-equilibrium systems. In the 
phenomenology proposed, no free energy can be weh defined. Fluctuations of Landau free energy 
couphngs arise when the intrinsic chemical potential leads to intrinsic disorder. The relaxation 
is shown to take the form of a geometrical random process. Systems of this kind show criticality 
features as well as that of first order transitions, which encapsulate in the form of a generalized 
static fluctuation dissipation relation. This will help determine three classes of distributions, which 
are, by defining a Hurst exponent for the relaxation rate: the regular Maxwell-Boltzmann for all 
if < |; the usual scale free universal distributions with power law tails for H G 1]; and a new 
class. The latter lies in the intermediate case, when H = ^. The distribution functions are scale 
free close to the origin, and cross-over to a Maxwell-Boltzmann asymptotic behavior, with both 
the scaling exponent and an effective temperature that depend on the magnitude of the intrinsic 
disorder. 
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I. INTRODUCTION 



Non-Gaussianity has often been reported, and has appeared universal in a broad series 
of experimental situations, ranging from turbulent flows to magnetic fluctuations, as, for 
instance, in the 2D XY model at low temperatures, and it has been related to criticality 

If one thinks of the Landau theory of fluctuations Maxwell-Boltzmann distributions 
describe monophasic fluctuations out of the vicinity of critical points, and for scales well 
above the correlation length. If metastable states are visited, it is in a way comparable 
to Arrhenius activation. On the contrary, in principle even at relatively low temperatures, 
out-of-equilibrium systems such as depicted in [ij] show a kind of disorder that displays swift 
and sharps changes with high magnitude. This means clearly that the relaxation process is 
not in a linear regime. 

Thus we need to seek out a different activation process. 

A. Random walks 

We start by following the remark that Maxwell-Boltzmann distributions can be taken 
from the limit central theorem by replacing time with the equilibrium time tgg (e.g. in a 
Ornstein-Uhlenbeck process) in the time dependent probabihty distribution t), 

p{6v; t) oc -y= exp • (1) 



It amounts to making a prescription for the collection of such particles 

{5v^) = t,,D, (2) 
which is the static fluctuation dissipation relation (FDT). This results in 

Peq[ov) OC -j== exp - . (3) 



From this, one already sees that it is not sufficient that tgq oo, but that regular 
diffusion ceases, in order to have non-Maxwellians. This means that the break down of the 
FDT relation eq.([2]) is necessary. 
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Another prototypical case, that one expects to happen close to criticality, is a scaling 
random walk, or Levy walk. Its time dependent distribution is of the following form[l] 

P{5v;t)^t-^g(^^'^ (4) 

with Q being a universal function. In just the same way, one takes the asymptotic limit 
t — 7- oo, and assumes that scaling extends to the collection. Then, normalizing with respect 
to time, one gets straightforwardly 

Pstat{Sv)<x5v-('+^). (5) 



So, there are at least two archetypes of probability distributions, depending on the nature 
of diffusivity in the system. One notices that the functional forms eq.([l]) and eq.(jl]) are 
identical, which is due to scaling between spatial journey versus time. However, there is a 
radical change in the physics, since, going from the one case H > ^ to the other < |, 
scale invariance ceases, and especially the correlation length becomes finite for a Brownian 
walk. 

The liquid gas critical point might be a good illustration. It is isolated in the phase 
diagram, and normal diffusion does not cease. One nevertheless observes opalescence, owing 
to the large fluctuations of the interface of the two phases. So, it lies in neither of the two 
already identified situations. 

The question that naturally comes to mind after this can thus be formulated as follows: 
is the functional limit with respect to H, of the probability distributions, in other words, of 
the universal function Q, continuous piecewisel The answer is no. 

In the following, I attempt to show what these kind of large fluctuations have to do with 
phase transitions of the first order, in interplay with the vicinity of a critical point, which 
can be of large extent. 

[1] Usually, one uses the fractional indexes a and /3, which equal respectively 1 and 2 in the Brownian case. 
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B. Slow dynamics 



On the one hand, break down of the FDT is a feature of slow dynamics. This term quahfies 
systems for which the correlation length depends on (usually grows with) time, generally a 
waiting time tm after some action, for instance, after the application of a magnetic field to 
a disordered magnet at low temperature j^. It may be written in the form 

it,,D,)=X(6^) (6) 



where X is a function, and where it is understood that the mean square depends slowly on 
time, 6v^{t^) Q, 5|. 

Now, if the large fluctuations are indeed associated with first order transitions, they imply 
latent heat exchanges, which means entropy exchanges. Thus, bearing in mind linear phe- 
nomenology, one may rightly suspect a global modification of dissipation, like a modification 
of FDT. 

On the other hand, large fluctuations of interest are stationary, as in 6|-l9|. Slow dynamics 
may become stationary if the external field is being randomly triggered, for example. If this 
is so, a random waiting time would then be associated with each random step, and the break 
down of the FDT would be averaged over the realization of such a process. It is likely though, 
from eq.(|6]), that a regular FDT would result from this operation, be it only from the central 
limit theorem. This, then, is where the influence of a critical point is essential, in providing 
a new specification to this random triggering: scale invariance (i.e. long correlations). The 
external field amplitude scales with the correlation length, and by means of their functional 
dependence for slow dynamics, it also scales with the waiting time. The only way that such 
a process can occur is if the external field is in fact intrinsic. 

We may now specify the form of the function X thanks to this hypothesis. 

II. AN EXAMPLE: BURGERS MODEL 
A. Extrapolation of linear phenomenology 

One dimensional systems are very unlikely in phase equilibrium, due to Landau's argu- 
ment [2!. Let us then consider a turbulent, compressible flow, in one dimension. There are 
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large fluctuations of velocity. Fluctuations of a certain magnitude have a certain spatial 
and temporal extension, therefore one can readily name phases. The situation mimics an 
inhomogeneous mixture of phases. With no other interaction, the chemical potential of each 
single phase is basically identifled in terms of the relative kinetic energy per particle 

/i = \5v\ (7) 

The thermodynamic current tends to equalize the potential of the different phases. Linear 
phenomenology relates it proportionally to the thermodynamic force 

F = -£V| 

with T = I {Sv"^), and C the phenomenological coefficient. The equation of phase motion is 
the inviscid Burgers mo del [2] 

^Sv = -X {V6v) 6v (8) 

with A = ^. 

As it is well known, a given perturbation may relax locally, but not globally towards some 
equilibrium. If a noise source term is added here, the system is incapable of equilibratin. 
with it, instead develops long range correlations [lo|, and no Gaussian distributions 



12|. 



So we can see that the most simple case falls out of the scope of the hypothesis underlying 
linear phenomenology and what it should be describing. The reason often alleged is that 
actually no thermodynamic potential can be defined for eq.([8]), so that the Burgers model 
has become a prototype for non-equilibrium dynamics. Let us now contest this fact. 



B. Dissipative function 

First, it is quite contradictory to treat a thermodynamic phase like a Langevin particle, 
by adding a noise source to the equilibration evolution of its chemical potential. In doing 
so, one mixes a microscopic modelization with a macroscopic description. A more coherent 
version should provide a macroscopic and self consistent mechanism that maintains the 

[2] With only one phase, = 0, the equihbrium of the fluid is given by Bernoulli's relation, which reads for 
a uniform density VT = ^S/Sv"^ 
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different 5v'^ phases out-of-equilibrium. It is noticeable also that the confusion might be a 
clue for seeking scale invariance linking micro to macro scales, in such a mechanism. 

Now, if one naively tries to find an effective free energy associated to eq.(|8]), it cancels 
out in the following way. If we look for a functional J-", such that j^J^ = A ( V5w) 5v, we 
try 

J^= j dx^^V5v^ (9) 
but then the variation is in fact zero, after integration by parts, 

=\ j dx^^ i5v) {\/6v'^ - \/5v'^) = 0. 

The cubic term cannot contribute to the variation of J-" because of its symmetries 
Sv — )■ —6v 

, besides having no absolute minimum. 

X — > —X 

We may nonetheless consider with good reason that the non local gradient and the local 
field are independent variables. 

Thus, instead of a cubic term we get an external field h for the intrinsic thermodynamic 
force in the free energy 

h = F = ^^Sv^ (10) 

There will be some link between the two conjugate fields, but only statistically. In addition, 
symmetries have to be preserved statistically, so 



(h) = 0. 



This in turn might yield (ST) = 0, and with no work exchange, at first sight, would seem 
to violate the second principle, but at this stage we should simply recall that the system is 
out-of-equilibrium. 

This treatment is a mean field description of the phases, because the intrinsic field h 
represents their large scale interactions that should lead to reorganization. It could be more 
general than the above expression eq. ffTOj) . Finally, the choice that was made is not a priori 
unique. Instead of h, a complex "mass" = could have been assigned 

ro = WSv 
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although it is less obvious then that its value adjusts self consistently, contrarily to what 
follows. 



III. GENERICITY OF FIRST ORDER FLUCTUATIONS: FLUCTUATING FREE 
ENERGY 

A. Degrees of freedom 

A first order transition close to a critical point resembles formally a subcritical pitchfork 



bifurcation, if one thinks of a dynamical system. Lorenz's famous system [13| is an example 
of the chaotic oscillations that occur in its vicinity. It could be that a particular phase, when 
forced from outside, be locally modeled by such a dynamical system, showing temporal chaos. 
But, it is not known what to expect when there is spatial coupling between different phases, 
because the number of degrees of freedom then scales with the system size. 

Schematically, just as it is not possible to identify a potential energy establishing a deter- 
ministic path for one single particle, because of the infinite number of microscopic degrees of 
freedom, but only to reduce to a statistical picture through a thermodjTiamic potential, the 
above infinite number of degrees of freedom, for an infinite system size, make unlikely the 
existence of a thermodynamic potential, which should have established a deterministic path 
for the collection of the particles. We point out here more precisely that the noise which is 
always used to model randomness and energy supply, in domain growth for instance, with 
Langevin models such as the Burgers model, should not be thermal noise, but pretty much a 
"non equilibrium noise", operating directly in the free energy rather than acting externally 
on the local order parameters or phases. 

In summary, an order parameter is indeed regarded as the state of a macro particle, 
under the influence of a random thermodynamic force. This means that no free energy min- 
imum truly exists, and that the macroscopic state diffuses out-of-equilibrium. The original 
Langevin model would have as a straightforward pending, something like a random set of 
thermostat baths in interaction with each other. 

The difference here is the necessarily self-generated nature of the noise. In a sense, one 
such macro particle must interact with itself, meaning that a macro particle designates a 
range from a single local phase to an encompassing of some phases. That is to say, a macro 
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particle understood as a set is in interaction with all of its subsets, and supersets. This is 
where scale invariance intuitively should come from. 



B. Renormalization orbits 

In Landau's theory, the external field h is considered to be a parameter, shifting the 
system away from absolute equilibrium, to an imposed equilibrium. If this parameter is 
intrinsic, as it is supposed to be here, it should enable the system to spontaneously leave 
absolute equilibrium. In this respect, it is valuable that the intrinsic field h comprises space 
in a hidden way, in conformity with the meta thermodynamic like approach we are looking 
for. 

Thus, within each set the value of h should be different, as if in a mixture of sets. This 
leads the values of the Landau free energy couplings to be different from set to set too. In 
order to make this practical, we use the ideas of scale invariance and renormalization 



16| . The way it works technically is simple: the free energy couplings undergo a subcritical 
pitchfork bifurcation, because of the finiteness of /i, and therefore the order parameter does 
so too. 

Consider initial conditions close to a critical state, and the effective Landau Free energy 
to be 3 

F = Feq + ^ro6v^ + B6v^ - hSv. 

The renormalization flow after rescaling by a factor b, for the parameters tq, B, and h, 
possesses one and a half degrees of freedom, as a dynamical system of equations. It is 
therefore likely to exhibit chaotic orbits as was noticed in [l6|. If so, because of stretching 
and refolding the orbit visits the neighborhood of the critical point as well as far off regions, 
besides it is likely that the orbit cycles chaotically from one side of the stable (fast ) ra ode 



variety to the other, therefore making the system go through first order transitions [17|. In 
other words, it shows hysteresis. It goes in a chaotic manner throughout each cycle from a 
critical state to a mean-field state. 

It should be noted that the argument is slightly different from the typical linear shell to 



shell decimation procedure 



18|, because it amounts to attributing to each single equili- 
brated phase a particular orbit cycle, with a typical cutoff correlation length, or extension 
length, associated with it. The cutoff (in correct units) corresponds to the inflexion point 
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from where the orbit returns to the vicinity of the stable mode variety. We then obtain for 
one single system sets of {^q}, {B^}, {h^}, and in particular {tg^}, where the superscript i 
stands for one of each cycle, indifferently over space and time. 

IV. A GENERALIZATION OF THE STATIC FLUCTUATION DISSIPATION RE- 
LATION 

The return point is far from the critical point, so mean field theory should be a good 
approximation there. At the same time, whenever the orbit goes by it, scale invariance 
should hold. So, a natural idea is that both phenomenologies hold together. 

A. Using mean field and scale invariance 

First, let us go back to the mean field theory. When tq = 0, the (non linear) susceptibility 
would be for h ^ 

= UBSv"^. (11) 
If we write instead of eq.([H]), loosing space dependence, 

dtSv = —h 

we would get with the help of eq. f lTT]) 

dt6v = -12XB6v^ 

but the solutions are not valid for all times. This is a reason to give to such a relation as 
eq. ffTT]) a conditional average meaning, with h being a random variable instead. 
Further, if we simply write the result of a diffusion process as 

{6v^) ~ D,t (12) 

where = J°°dT{hhr), because we would like a stationary process, it is essential that 
anomalous dimension arises in order to cancel out the secular divergence. In particular Dyt 
should be taken as a whole, and depend on the value of 6v'^ in the form of eq.([n]). 
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Now that we have a mechanism which is self generating noise, we prefer to model the 
relaxation process with a stochastic 'mass' = vq, that we write (£ is taken constant with 
no loss of generality) 

Aro = A. (13) 
So relaxation in the system is modeled by a geometric process 

dtSv = -~X6v. (14) 

For instance, if A had no relation with 6v, a finite variance, and was short range correlated, 
eg ■ 0141) would be a multiplicative cascade, leading to a non stationary log-normal distribution 



19|. Define a Hurst exponent coming from the conditional correlation function of A as 

dT t dT (\{t)\{t)/5v'^ ^t'''. 



By integration of eq. (fT4|) we obtain 



2H 



but a two time relation, reminding us of the slow dynamics, can also be found at a time r, 
by integration over a much smaller time interval t 



{6v\t)) (t) 6v^{T,t)t^'' . (15) 

This can be regarded as intermediate thermalized states reached at times r. We obviously 
want a real constant {Sv"^), defining the amplitude of the stationary fluctuations. As already 
mentioned, the r.h.s. of eg. f|T5|) is therefore considered inseparably, meaning that (dropping 
over line notation) 

Whereas ro — )■ 0, this asymptotic in time prevents the averaged susceptibility to become 
infinite (Tx) oo. We then can work it out by replacing outrightly t^^ with A, very much 
as we have done so, in an eguilibrium state, for the correlation time in eg. ([2]). This is almost 
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by the definition of tfie set of all possibilities for this random variable 

Therefore, the fluctuating relaxation rate scales as 

'\/5v)^6vTi. (16) 



In particular {D^t) is now a constant, and the function X[A) oc A""^ . 

Actually, for the phenomenology of orbit cycles to hold, one needs a non zero overall 
dissipation (or instability) Aq = Aro 7^ 0, although Aq — )■ 0. Thus, eq.f lT6|) is rather 

(\/5ii^ ^ Ao ± C±5vTi (17) 

where C± is related to B. For instance, for a conserved quantity it describes locally the 
effects of the energy transfer from scale to scale within a physical range, where the effective 



diffusion rate is -D,, oc A^ 



2Q| 



Finally, it is worth noticing that eq. (|T7|) renders a kind of conditionally averaged Landau's 
free energy, of the original form as soon a.s H = i, which is indeed the case here. Otherwise 
because of eq.f ll3p . the coherence of the phenomenology starting with Landau's free energy 
might be questioned. 

B. A guesstimate on indexes 

Because near the critical point, it is the linear behavior that dominates, we can assume 
that 



X = h^-'g — (18) 

where tq oc e, and g{x) — )■ as a: — > 00. The interesting limit here is a: — )■ 0, where to first 
order 

\img[x) = cte. 
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Now, by identification of tlie scale factors in eq. fllSI) and eq. fll6l) . one finds tlie Hurst exponent 

H=^. (19) 

7 

If the phenomenological coefficient C is not anomalous (what we have supposed here) then 
we also have 

7 = zu. 

One has to be careful though that the critical indexes may not have their regular values. 
Indeed, in the linear hyperbolic region one can write tqA' ~ 1, where A stands for the 
eigenvalue and / for the iteration integer, as soon as the correlation length scale is reached. 
For one particular cycle of the chaotic orbit, the iteration integer starts from zero, and one 
can rewrite this as 

r^A'» = a\ 

Letting down the superscripts for clarity, noting the correlation length ( = UQ, where Q is 
a fixed number, if it happens that 

a ~r[^ 

then the critical index v will be modified by a factor (1 — a) as 



Incidentally, it has been suggested in [17| that the broken symmetries due to h as well as 
the related disorder could change the indexes' values, for instance from a 'bosonic' to a 



'fermionic' mean field, for which v = 1 and 2 = 1, as was observed in 

Leaving this indetermination aside for the moment, interestingly, when H = ^ the system 
seems to be at the border of normal fluctuations, and from eq. lfTOl) . it seems to be related 
to the mean field values of the indexes (3 = ^, u = ^ and z = 2 (or u = 1 and z = 1). So, 
H = ^ could be related to marginality. 
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V. STATIONARY DISTRIBUTIONS 



The previous generalization of the FDT allows the derivation of explicit forms for the 
probability distributions. 



A. Functional form 



Consider first when G ] |, l] , we have for the Fokker Planck criterion 



^2H 

lim = 0. 

t^o t 



Thus the Fokker Planck equation is reduced to = — ^ (^(^\/6v^ 6v . This leads 
to the following stationary distribution, p(Sv) oc exp (— J dSv j) , where the integrand is 

(\/Sv\ -\-Sv-Jy- (X/Sv\ . I . r> 1 •! r 1 T -1 • • 

I = ' /7,7\ ■ Dj using eg. ( Hop one nnds a power law tail for the distribution m 

dv(\/ov) 



accordance with eq.(l2T|) 

p{6v) (X 6v^(^-^f\ (21) 



The exponent 



zu 



is called the Gutenberg Richter exponent, and belongs here to the interval ]2, 3]. 
If < |, one gets the usual Maxwellian distribution eq.(|3]). 
Similarly, the intermediate case if = | results in 

Pstat {Sv, (P)) OC 6v-'^^ exp I^^C^-^6v'^ . (22) 

Close to the origin, the behavior is dominated by a power law. Interestingly, the slopes rely 
not only on the combination of the indexes 1 + jj, but also on the noise intensity ^A^ 
Asymptotically, one recovers a Maxwellian tail, with an effective temperature 



Teff^[X'). (23) 
The cross-over goes to infinity with the effective temperature. Here, we have put back an 
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overall dissipative 'mass' Aq, and the generalized relation, accounting for signs, looks like 
eq.([I7D 



VI. CONCLUSION 

I have presented a phenomenology for the determination of the statistics of out-of- 
equilibrium fluctuations, based upon the existence of a free energy, but where its couplings 
cannot be defined properly. These are indeed likely to cycle chaotically. Systems so governed 
go through second and first order transitions at the same time. A distribution (which resem- 
bles a chi square only in a particular case) that reminds us of both a scale free distribution 
at small amplitudes and a Maxwellian at larger amplitudes, with an effective temperature 
expressing the strength of intrinsic disorder, is found. 

The interplay between first and second order transitions was first found in the study of 
the hysteretic behavior of some fluids, where it is due to instabilities of the Rayleigh type. 
I only mention here a common fluid modelization for avalanches in sand piles and for the 
resonance between waves and particles in a plasma [l7|, as this will be considered elsewhere. 

I now would like to evoke further work that I think is worth while investigating, given 
this phenomenology. For instance, J. Sethna et al. [22| have also found a cross over behavior 



for the PDFs of avalanche amplitude in the random field Ising model, which bears a priori 



extrinsic disorder. S. T. Bramwell et al. 



have performed precise calculations for the 



spin wave approximation of the XY model which establish evidence of a relation between the 
different moments for the Fokker Planck perturbative development. This could be related 
to some kind of Noether (or Ward) identity, in marginal systems, due to scale invariance 
symmetry, that we have set out here as a modified static FDT. I also question whether the 
phenomenology of chaotic cycles may help in the understanding of intermittency, which lies 
in the UV range, in developed turbulence. 
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